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Performance Statistics of the FORTRAN IV (H) Library 
for the IBM System/360 


A report describing testing procedures and results 
for accuracy and timing tests of the basic IBM 360/50 
FORTRAN IV (H) subroutine library has been pre- 
sented. The testing was undertaken to verify perform- 
ance capability and as a prelude to providing some 
replacement routines of improved performance. 

All certification runs were made under the standard 
operating system, OS/360, on a Model 50 computer. 
The tests were divided into three categories: accuracy, 
timing and error returns. 

The greatest effort was directed toward the ac- 
curacy tests. Consider the testing of a single-precision 
function subroutine. First the domain of the routine is 
divided into several subdomains; the subdomains usu- 
ally correspond to subdomains used internally in the 
routine under test. Next, several thousand uniformly 
distributed pseudo-random, single-precision, floating- 
point numbers, together with their identically valued 
double-precision counter-parts, are generated for each 
subdomain. The single precision numbers are used as 
arguments for the routine being tested; the double- 
precision numbers are used as arguments for the cor- 
responding double precision subroutine. The results of 
the double-precision computation are considered as 
“correct.” The single-precision results are then con- 
verted to double precision and the maximum relative 
error and the root-mean-square of the relative error 
are computed in double precision. In addition, the 
single-precision results are compared against the re- 
sults of the double-precision computation correctly 
rounded to single precision, and the frequency of the 
errors in units of the last bit of the mantissa is tabu- 
lated. Similar tests were performed on the double- 
precision subroutines. 


Timing tests were performed by invoking the rou- 
tine under test up to 10,000 times, with either a fixed 
argument or random arguments, in a loop. The “over- 
head” of the looping instructions was timed in each 
case by repeating the same loop with a dummy sub- 
routine whose only executable instruction was a jump 
back to the calling program. Because of the coarseness 
of the internal clock, the indicated timings varied by 
as much as 5 m sec. from run to run on a given test. 

The tests of error returns consisted in supplying 
arguments outside the domain of the subroutine under 
test. Without exception, this resulted in termination 
of the job. 

The report includes tabular summaries of the above 
tests, a listing of the decimal equivalents of selected 
inverse powers of 2, and a typical test program. 
Notes: 

1. The report is entitled “Performance Statistics of 
the FORTRAN IV (H) Library for the IBM 
System/360,” ANL-7321, May 1967. This report 
is available from the Clearinghouse for Federal 
Scientific and Technical Information, Springfield, 
Virginia 22151. The printed copy is $3.00 (micro- 
fiche copies, $0.65). Additional details will be 
made available in “Performance Testing of Func- 
tion Subroutines” by W. J. Cody, Jr., Proc. 1969, 
S.J.C.C. Paper to be presented at Spring Joint 
Computer Conference, Boston, May 16, 1969. 
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Patent status: 

Inquiries concerning rights for commercial use of 
this innovation may be made to: 

Mr. George H. Lee, Chief 
Chicago Patent Group 
U.S. Atomic Energy Commission 
Chicago Operations Office . 

9800 South Cass Avenue 
Argonne, Illinois 60439 
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